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ABSTRACT 

We argue that massive stars are the dominant sources of energy for the turbulent motions 
within giant molecular clouds, and that the primary agent of feedback is the expansion of Hll 
regions within the cloud volume. This conclusion is suggested by the low efficiency of star 
formation and corroborated by dynamical models of Hll regions. We evaluate the turbulent 
energy input rate in clouds more massive than 3.7 x 10^ solar masses, for which gravity does 
not significantly affect the expansion of Hll regions. Such clouds achieve a balance between the 
decay of turbulent energy and its regeneration in Hll regions; summed over clouds, the implied 
ionizing luminosity and star formation rate are roughly consistent with the Galactic total. Hll 
regions also photoevaporate their clouds: we derive cloud destruction times somewhat shorter 
than those estimated by Williams and McKee. The upper mass limit for molecular clouds in 
the Milky Way may derive from the fact that larger clouds would destroy themselves in less 
than one crossing time. The conditions within starburst galaxies do not permit giant molecular 
clouds to be supported or destroyed by Hll regions. This should lead to rapid cloud collapse 
and the efficient formation of massive star clusters, which may explain some aspects of the 
starburst phenomenon. 

Subject headings: ISM: clouds - H II regions - stars: formation 



1. INTRODUCTION 

Giant molecular clouds (GMCs) are the sites of 
most star formation in the Milky Way, and the evolu- 
tion of the Galaxy and its stellar population are con- 
trolled in large part by the physics of GMCs. Despite 
several decades of observations, the clouds' physical 
nature, longevity, and modes of formation and de- 
struction are still matters of debate. If clouds survive 
for more than a single dynamical time, it is plausi- 



ble that their observed properties are (like stars') the 



product of internal sources and sinks of energy (for a 
recent review, see McKee 199E). If instead the y dis- 
appear after one crossing time ( Elmegreen 2000|) then 
their properties will be more affected by the mech- 
anism of their formation (|Vazquez-Semadeni et al 



1.1. Cloud Scaling Laws and Their Origins 

Molecular clouds are observed to obey a set of scal- 
ing relations known collectively as Larson's ( |1981 ) 
laws. As updated by ^olomon et al. (1987 ), these 
can be summarized: a constant mean column density 
S ~ 170AfQ/pc^, corresponding to a visual extinction 
Ay ~ 7.5 mag; and virial balance, with a virial pa- 
rameter a = 10 ^rms/(3^esr.) (|Bertoldi fc McKee 1992|) 
of order unity ( [Myers &: Goodman 1988 ). Solomon 
et al. adopt the value a = 1.11, yielding the relation 
Urms — Vesc/V^ between the r.m.s. and escape veloc- 



1997; Ballesteros-Paredes et al. 1999a, b^^ However 



ities. A third but not independent property is the 
scaling of line width (cj) with cloud radius (i?ci = -Rpc 

/ — 1 /2 

pc): o" = i^rms/v3 = 0.72i?pc km/s. As they do not 
represent infall, these motions are considered turbu- 
lence. The molecular gas is cold (~ 10 — 30 K) and 
highly magnetized (~ 30/xG), so its motions are su- 
personic but roughly Alfvenic. 

Vazquez-Semadeni et al. (1997) advance a model in 
which these properties derive from motions within the 
interstellar medium (ISM) from which GMCs form. 
In two-dimensional simulations of turbulence within 
the Galactic disk, they identify a population of over- 

^We only consider clouds like those in the ^olomon et al, survey, and assume this survey correctly infers the clouds' properties. 
As we shall concentrate on the most massive clouds (Aici > 3.7 x 10^ M©), whose observations may be affected by beam smearing 
and velocity crowding (J. P. Williams 2001, private communication), this assumption should be checked by future observations. 
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rapid destruction also implies that the means of cloud 
destruction may play an equally important role. This 
paper considers what consequences the mechanisms 
of cloud destruction - primarily Hll regions - imply 
for the internal motions, energy budgets, star forma- 
tion rates, and lifetimes of GMCs. 
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densities that could be considered clouds. Because 



of the turbulent spectrum, all of the objects formed 



who implicated main-sequence winds, and Lada & 
Gautier (1982), who realized the potential of their 



in their simulations obey the line width-size rela- 
tion noted above. However, all but the few most 
massive of these are transient compressions rather 
than self-gravitating objects; hence their escape ve- 
locities and surface densities do not follow the virial 



relati pn. [Vazquez-Semadeni et al.| suggest that se^ 



lection effects restrict observed clouds to a narrow 
range of inferred column densities. Specifically, they 
note that the IRAS survey of [Wood et al. (1994| ) 



may only be sensitive to an outer shell of warm dust 
around clouds of various columns. However, |McKeg| 
(1999) has countered that CO would form and be de- 
tectable at significantly smaller columns {Ay = 1.8; 
[van Dishoeck &: Black 1988| ) than are typical of the 
observed GMCs. 

An alternative possibility is that the common CMC 
column density arises from internal cloud processes: 
specifically energetic feedback due to star formation 
(as in the model of [McKee 1989] , hereafter M89). For 
this to be possible, two conditions must hold: 1. stars 
form rapidly only in regions that exceed a critical col- 
umn density; and 2. star formation is potentially so 
vigorous a source of turbulent energy that it over- 
whelms the natural decay of turbulence if the column 
is much above this critical value. Under these condi- 
tions a GMC will settle into a state of energetic equi- 
librium much like a star's, with star formation occur- 
ring just fast enough to offset turbulent decay. The 
necessary column density would be roughly the crit- 
ical value, hence the common value observed among 
GMCs. 

For condition (1), [McKee proposed that star for- 
mation is inhibited in regions that are not shielded 
by ~ 4 visual magnitudes in each direction, because 
these layers are penetrated by far-ultraviolet (FUV) 
photons that elevate the level of ionization. FUV ion- 
ization thus slows ambipolar diffusion, which M89 ar- 
gued to be a rate-limiting step for star formation be- 
cause stellar mass regions are too highly magnetized 
to collapse directly. This sets the critical column den- 
sity at roughly 8 visual magnitudes (4 on each side), 
close to the value of 7.5 observed in GMCs. Corrob- 
orating this hypothesis, Onishi et al. (1998| ) find a 
sharp distinction among substructures in the Taurus 
clouds between those with Ay ^ 7.5 mag that are ac- 
tively forming stars, and those with lower extinctions 
that are not. In the low-metallicity environment of 
the SMC, pak et al [T998[) verified that the clouds 
maintain Ay ~ 7.5 mag, although this required a 
higher column density than for Milky Way clouds. 

For condition (2), McKee[ specified protostellar 
winds as the agents of star formation feedback, follow- 
ing |Norman^SilF(198^) and [Franco k Cox (1983[) , 



protostellar counterparts. Subsequent numerical sim^ 
ulations (e.g., Mac Low et al. 1998[; ptone et al. 1998 ; 



Mac Low 1999[ ; pstriker et al. 2001[ ) have indicated 
a much faster decay of turbulence than McKee as- 
sumed, calling into question (e.g., Basu fc Mural i 
200 1[ ) the notion that energy injection from protostel 



lar winds is vigorous enough to offset turbulent decay. 
Below, we show that Hll regions (originally consid- 
ered by Arons Sz Max 1975[ ) represent an additional, 
inevitable, and more important source of turbulent 
energy for molecular clouds. 

1.2. Molecular Cloud Destruction and the 
Inefficiency of Star Formation 

The rapid decay of turbulent energy should logi- 
cally cause an equally rapid contraction of GMCs, the 
end product of which is star formation. Nevertheless, 
clouds make stars at a tiny fraction of the rate al- 
lowed by dir ect gravitational collapse ( pucker man & 



Palmer 1974 ). Star formation is not only slow, but in- 



efficient: the protostellar sources observed by Cohen 
fc Kuhi (1979|) compose ~ 10% of their surrounding 



clouds; Myers et al. (198^ ) surveyed 54 clouds, esti- 
mating ~ 2% of their mass to be in stars, on the basis 
of their Hll regions; and [Wilhams fc McKee (1997[ ) 
(hereafter, WM97) argued that only 10% of the mass 
of a GMC would ever become stellar. More recently, 
[Carpenter (2000D has estimated that 1% - 9% of the 
mass of several nearby molecular clouds is in embed- 
ded stars. The sluggishness and inefficiency of star 
formation can only be consistent with the rapid de- 
cay of turbulent energy if GMCs are destroyed more 
rapidly than they can convert themselves into stars. 

What process destroys m olecular clouds? 
[Ballesteros-Paredes et al. (1999b|) note that gravi- 
tationally unbound clouds formed by turbulent com- 



pressions (most of the objects seen by [Vazquez- 



Semadeni et al. 1997) are easily disrupted in a single 



crossing time by the flows that created them. How- 
ever they also show that bound objects, once formed, 
continue to collapse rather than re-expanding. (In 
these authors' simulations, such clouds are stabilized 
or disrupted by a local heating included to represent 
the action of massive star formation, a topic to be 
addressed in detail below.) The Solomon et al. sur- 
vey indicates that molecular clouds are far too tightly 
bound to be destroyed by turbulence in the interstel- 
lar medium: their hydrostatic pressures are in excess 
of 3 X W^kB cm , whereas the bounding gas pres- 
sure is only 2 x lO^ks cm~'^ ( [McKee 1999| ). The ram 
pressure of motions in the diffuse ISM (uh ~ 1 cm~^, 
fj ~ 10 km s~^) is also ~ 2 x lO'^A;^ cm~^, and thus 
insufficient to disrupt a cloud. Clouds contain a bind- 
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ing energy per volume that is roughly three halves 
their hydrostatic pressure (somewhat less, due to 
magnetization; McKee 19991 ); this was radiated in 
the process of formation and must be resupplied to 
unbind the cloud. The most plausible source for this 
energy is young and massive stars formed within the 
cloud itself. 

Several mechanisms have been considered in this 
regard. Matzner Sz McKee (2000| ) have shown that 
small clouds can be disrupted by protostellar out- 
flows as they form low-mass stars; however, they pre- 
dict such mass loss not to be important for giant 
clouds. Elmegreen (1983| ) suggested that radiation 
pressure might unbind clouds if their stellar popu- 
lations became too luminous, but neglected the re- 
processing of radiation into far-infrared wavelengths 



to which clouds are transparent (see Jijina fc Adam: 
1996|]|. Molecules are destroyed by FUV photons in 
photodissociation regions; however, the thermal ve- 
locities of these regions are far below the escape veloc- 
ities of GMCs. Therefore, photodissociation regions 
are thought to form an atomic layer around GMCs 
( |van Dishoeck &: Black 1988 ) and have only been sug- 
gested as a means of destruction for clouds that are 



dynamically disrupted by a different process (Whit 



wort lj 1979|; [McKee &: Wihiams 1997|; [Hartmann et al 



2001|)|, as lWilhams &: Maddalena (1996|) " may have ob- 
served around the cold cloud G216-2.5. 

Whitworth (1979) calculated the destructive effects 
of massive stars due to the ejection of photoionized 
gas, finding ~ 10^ Mq to be ejected for each blister- 
type Hll region. Supernova explosions typically add 
~ 10% and at most 40% to this amount ( [Yorke et al. 
19891) , for HII regions created by single stars. A 



few supernovae arise from B stars in the mass range 
8 — 12 Mq, which do not have appreciable Hll re- 
gions ( [Chevalier 1999| ); however, these require at least 
1.3 X 10' years to evolve, and eject only ^ 500 Mq 
each if they explode inside their cloud (and much 
less otherwise; Tenorio-Tagle et al. 1985]) . The main- 
sequence winds of massive stars emit a total energy 
comparable to supernovae. McCray fc Kafatos (1987 ) 
and McKee et al. (1984 ) have argued that stellar wind 
bubbles are confined within Hll regions in the con- 
text of a cloudy medium. For these reasons, we con- 
centrate on photoionization as the primary cause of 
destruction for GMCs (as have Whitworth 197£ 



Bhtz 



|fc Shu 198q , [Yorke et al. 1989| , [Franco et al. iggif aJid 
WM97). The dynamical effects of stellar winds and 



supernovae are considered in §3.2 



The fact that massive stars are responsible for un- 
binding molecular clouds more rapidly than they can 
form stars implies that they supply more energy than 
is dissipated in turbulence. If a fraction of this en- 
ergy is incorporated into the motions of the remaining 



molecular gas, it will sustain turbulence and slow the 
cloud's contraction. In the following sections we con- 
sider the loss and regeneration of turbulent motions 
(^, the relative importance of Hll regions and pro- 
tostellar winds (§ [2.1[ ), the dynamics of an individual 
association (§^) and of a population of associations 
(^, and the implications of an equilibrium between 
the decay of turbulence and its regeneration in Hll 
regions for the ionizing luminosities and star forma- 
tion rates (^ and lifetimes (§^) of massive clouds. 
These estimates are corrected for the interaction be- 
tween Hll regions in 

Finally, in §^ we point out that the high-pressure 
environments of starburst galaxies prevent molecu- 
lar clouds from being destroyed or even supported 
by photoionization. In the absence of other sources 
of energy, such clouds must collapse and form stars 
at a high efficiency; this is a recipe for the efficient 
formation of massive star clusters often observed in 
starbursts. 

2. SOURCES AND SINKS OF TURBULENT MOTION 
To consider the fates of GMCs, one must account 



for the gains and losses of turbulent energy ( Norman 



fc Silk 1980| ) . M89 assumed a range of timescales for 
the dissipation of turbulent energy. 



Etmh 



turb,diss 



(1) 



ranging from three to ten times the free-fall time of 
the cloud. Here, -Eturb represents both kinetic and 
magnetic energy associated with turbulent motions 
([Zweibel &: McKee 1995 ) . Simulations by [Stone et al 



|(1998| ) show incomplete equipartition between these 
components: for the case that most resembles molec- 
ular clouds (magnetic pressure 100 times greater than 
gas pressure). 



turb 



0.78Mdt;2^3. 



(2) 



McKee s assumption of a relatively long decay time 
scale was based on the notion that the magnetic field 
should cushion gas motions. Testing this notion nu- 
merically, several groups ( Mac Low et al. 1998 ; Stone 
et al. 19981 ; [Mac Low 1999[ ; [Ostriker et al. 2001|) have 
found the dissipation to be much more rapid, espe- 
cially if turbulence is driven on scales smaller than 
the entire cloud. For the same physical situation that 
gives equation (Q), [Stone et"ar find 



0.83- 



(3) 



where Ain is the wavelength on which the turbu- 
lence is stirred. Combining this with [Solomon et al 
(1987[ )'s virial relation v^ras = fesc/\/3 gives tdiss/^flf — 
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0.93Ain/-Rci- This dissipation time scale is shorter 
than the range assumed by M89 if the forcing scale 
for turbulence is smaller than the radius of the cloud; 
the importance of the forcing scale has recently been 
highUghted by |Basu &: Murah (20011 ). 
To estimate how turbulence is driven one must al- 



low for the radiative nature of the gas, which causes 



compressions to be very dissipative. Our argument 
folloA ^s [Norman fc Silk (198q) and M89. An impulse 
(such as a protostellar wind) whose momentum is 6p 
will cause a disturbance that decelerates as it sweeps 
into the cloud. Since energy is radiated, a thin shell 
forms and conserves linear momentum in each direc- 
tion (e.g., Matzner Sz McKee 1999a] ); the kinetic en- 
ergy associated with the motion is v5p/2 when the 
velocity is v. This continues until v has decelerated 
to a terminal velocity that WM97 estimate to be the 
effective sound speed, Cgg = \fPjp. (Note that Cgs 
far exceeds the thermal sound speed, as P includes 
the total hydrostatic pressure.) At this point the 
swept-up shell thickens, stalls, and loses coherence, 
rendering its energy to the turbulence. The increase 
of turbulent energy is therefore 



(4) 



The efficiency coefficient is uncertain, and must 



be determined by simulation; [McKee (1999]) suggests 
^ ~ 1.6 to account for the energy stored in magnetic 
perturbations at the end of deceleration. Equation 
(Q) implies that protostellar winds and Hll regions 
generate turbulent energy in proportion to the rate 
at which they impart momentum to the cloud. 

Equation (^) is ess entially the same form ula em- 
ployed previously by Norman fc Silk (1980|), M89, 



Bertoldi h McKee (1996|), pVIcKee (1999| ), [Matzner & 
McK jie (1999b[) and [Matzner (1999|) . It assumes that 



energy is injected in an explosive manner, so that the 
early stage of momentum-conserving thin shell ex- 
pansion (at speeds above Cgff) can be separated from 
the later stage of turbulent dissipation (at speeds of 
about frms)- It is only valid in cases where the kinetic 
energy of relative gas motions is present in the center- 
of-mass frame of the cloud; therefore, it does not ap- 
ply to large-scale gravitational fields. For the same 
reason its validity is suspect if the force that gives rise 
to S-p is applied over a large length scale compared to 
the cloud radius. However, we shall see in the subse- 
quent sections that Hll regions are explosive events 
that input momentum on scales smaller than i?ci, at 
least for clouds in our mass range of interest. 

What are appropriate values of bpl In the forma- 
tion of a low-mass star, a fraction of the material 
that accretes onto the star-disk system is redirected 
into winds with a characteristic velocity ~ 200 km/s. 



An observational analysis by Richer et al. (2000 ) 



implies a characteristic wind momentum of roughly 
50 km s~^ times the mass of the star that forms; how- 
ever, this analysis is quite uncertain. Models exist in 
which the wind removes anywhere from a tenth to 
half the mass flowing through the disk ( Pelletier & 
Pudritz 1992[ ; [Najita fc Shu 1994[ ). We adopt one- 



sixth as an intermediate value (following [Matzner & 
McKee 2000| ), so that the wind mass is five times 



smaller than the star's mass. With this choice, each 
star of mass generates a wind impulse 



,,171^, X 40 km s 



(5) 



where (/>^ represents our uncertainty. In § [3.2| we find 
that main-sequence and evolved stars contribute a 
comparable momentum in the six million years after 
their formation; here, we restrict our attention to the 
more impulsive protostellar winds. 

If is the total rate at which mass is converted 
into stars in a cloud, equation (H) gives 



turb.ui 



-M^(40 km s-^)ceff 



(6) 



as the rate of turbulence regeneration by protostellar 
winds. 

Next, consider cloud destruction at a rate Mdest 
in a sequence of discrete events (ejecting SM^est per 
event), no one of which completely destroys the cloud. 
Each mass ejection delivers an equal and opposite 
impulse to the cloud material left behind, and each 
impulse increases cloud turbulence according to equa- 
tion (m. S ince the ejection velocity must exceed fesc 
near the surface of the cloud. 



(7) 



If photoionization is the primary means of destruc- 
tion, then material is ejected at about the thermal 
velocity of ionized gas, cn ~ 10 km s~^, and 

Spdest = 4>I1 (^MdestClI, (8) 
where (pu is defined in analogy to (/>^; by equation 



turb,dest 



- -/Widest CllCeff- 



(9) 



An analysis of the generation of momentum in ion- 



ization fronts, presented in section [3.1[ , shows that 
011 > 2. 

2.1. Massive Stars Dominate the Energy Budgets of 

GMCs 

A comparison between equations (^) and @ re- 
veals the importance of massive star formation in the 
support of GMCs. Relevant to this comparison is the 
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net efficiency of star formation, SFEtot ) defined as the 
fraction of cloud mass that wih ever become stars: 



SFE 



tot 



/ M^dt 



(10) 



the denominator equals the initial mass of the cloud 
(in the absence of ongoing accretion from the ISM). 
Thus 



est 



(SFE 



tot 



-1 



DM. 



(11) 



where the brackets are time averages. SFEtot is esti- 
mated observationally by the current fraction of mass 
in stars; the two are comparable if the efficiency is 
low, if stars remain within their clouds ( Matzner &| 
McK jie 2000]) . If we take SFEtot ~ 5% as typical, then 
95% of the cloud's mass will be disrupted by blister 
Hll regions rather than achieving stardom. Equa- 
tions (^), (^, and (11) imply that the momentum 
contributed by protostellar winds is smaller than the 
contribution by Hll regions so long as 



SEE 



tot 



< 



1 



4(/'^/0ii + 1 



(12) 



Clouds that make stars as inefficiently as Galactic 
GMCs derive more turbulent energy from blister Hll 
regions alone than from protostellar winds. In §Q we 
shall show that HlI regions would be more important 
than protostellar winds even if star formation were 
efficient. 

3. EFFECT OF A SINGLE ASSOCIATION 

3.1. Momentum Generation by an Hll Region 

The dynamical phases of expansion of an Hll re- 
gion have been presented in numerous prior works 
(e.g., Whitworth 1979| ; Spitzer 1978| ). These authors 
neglected the inertia of the shell of shocked cloud 
gas; we give only a cursory treatment including this 
inertia, for the purpose of quantifying the momen- 



tum generated. We follow McKee &: Williams (1997 ) 
(hereafter, MW97) in adopting a common tempera- 
ture of 7000 K for ionized gas (and mean molecular 
weight 0.61, so cn = 9.74 km/s) and in approximat- 
ing that 27% of the ionizing photons (emitted at a 
rate IO'^^S'49 per second) are absorbed by dust rather 
than gas. We will also take for the molecular gas 
a uniform density po and hydrogen number density 
uh consistent with a mean hydrogen column density 
Nh,22 = iVj^/(1022cm-2) = 1.5 (ISolomon et al. 1987|) . 

After a rapid expansion to the initial Stromgen ra- 
dius (using the recombination coefficient of ^torey ^ 
Hummer 1995| ) 



Rst,o = 2.9 



,22 



1.5 



" sII'mX! PC, 



(13) 



the Hll region is governed by the requirement that 
a very small fraction of the ionizing photons actually 
reach the ionization front. This causes ionized gas 
density pu to vary as 



Pu 
Po 



Rst,o 



-3/2 



(14) 



where rn is the current radius of the ionization front. 
As the ionization front expands subsonically the den- 
sity Pu and pressure piiCjj are nearly uniform within 
the HlI region (albeit more perfectly for embedded 
than for blister regions, in which a pressure gradient 
develops as gas accelerates away). 

So long as the Hll region expands supersonically 
with respect to the molecular gas, it is bounded by 
a thin shell of dense, shocked cloud gas. The ra- 
dius of this shell is nearly identical to rn, and we 
shall restrict our attention to the period of expansion 
ni ^ Rstfl, when nearly all of the mass originally 
within rn remains within the shell (Spitzer 1978). 
Let Ash denote the shell's area and Mgh its mass. 
The shell's momentum equation is 



di 



(Mshni) = AshPii Cn + nii(nii - hi 



(15) 



On the left of this equation is d5p/dt, the rate of 
increase of the shell's momentum; on the right, the 
forces due to pressure (first term in brackets) and due 
to thrust caused by the exhaust of ionized gas at a 
velocity uu relative to the cloud. 

Blister and embedded Hll regions differ in the co- 
efficient relating A^h to rn: ^sh = (1)2) x 27rrjj 
for (blister, embedded) HII regions whose ionization 
fronts are idealized as hemispheres and spheres, re- 
spectively. 

Another difference is the relative importance of 
thrust and pressure in generating momentum. In 
an embedded region, the ionized gas is trapped and 



only expands homologously; equation (14) then im- 
plies Uu = rn/2, which is much less than cn when 
rn ^ Rst,o (iSpitzer 1978 ) so that only pressure need 
be considered. In blister regions, on the other hand, 
ionized gas flows away freely allowing the ionization 



front to tend toward the D-critical state (Kahn 1954) 
for which uu — hi = —cu if the ionized gas is effec- 
tively isothermal: recoil is just as important as pres- 
sure in generating momentum. The term in brack- 
ets on the right hand side of (15) is therefore, to a 
good approximation, (2, 1) x Cjj for (blister, embed- 
ded) regions respectively, assuming a D-critical ion- 
ization front for the former. 

The rate at which mass is ionized is puihi—uu)Ash, 
or approximately pn^shl""!!! for a blister region when 
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ni ^ Rst,o- The ratio between the rate of momen- 
tum generation and the rate of ionization is (pucu, so 
equation ( [T5| ) imphes 



cii 
-nil 



+ 



cii 



> 2 



(16) 



-Cii, i-e., a D-critical 



where equahty holds for uu 
front. 

The expansion of the Hll region is simplest to de- 
termine when rii ^ -Rst.o- In this phase, Mgh — 
^shPoni/3 assuming radial expansion. We seek a self- 



similar expansion of the form rn oc f; equation (15) 
admits the solution^ 



(2,1) X 3 
rj{4rj — 1) 

(23, 19) X 



2/7 



3.7 Myrj 



N 



H, 22 

I. 5 



-3/7 



xM^/^9^/^ DC 



(17) 



for (blister, embedded) regions, respectively. In this 
equation the first line indicates r] = 4/7, which fixes 
the coefficient in the second line. We have normal- 
ized to the ionization-averaged lifetime of rich OB 
associations (M97) of 3.7 Myr. 

The momentum of radial motion of the expanding 
shell is 



6p 



(2.4,2.2) X 10' 



t 



3.7 Myr/ V 1-5 



xM^/^^sti'^ Mq km s-^ 



(18) 



for (blister, embedded) regions. Note that the ex- 
tra thrust generated at the ionization front in a blis- 
ter region compensates for the smaller working sur- 
face Ash- Indeed, the two results are so similar that 
we may estimate 5p using the intermediate coefficient 
2.3 X 10^ Mq km s""*^, without discriminating between 
the two types of regions. This will simplify the anal- 
ysis in §^. 

For a blister region, the mass evaporated is 
5p/{(j)iicu): for cpu = 2, 



SMdest = 1.2 X 10^ 



t 



3.7 Myr 



''H,22 

1.5 



(19) 



favorable comparison gives us confidence in the ap- 
proximations adopted in this section. 

Note that the above results differ quantitatively 
from those given by [Whitworth (1979| ), which were 
adopted by WM97 to study the erosion of molecu- 
lar clouds (see §3T). Our equation (|l9|) for the mass 
evaporated agrees with these authors' results within 
1%. However, equation ([TtI ) gives a radius for a blis- 
ter or an embedded region that is smaller by a factor 
of 1.6 or 1.9, respectively, than the characteristic size 
quoted by these authors for blister regions - imply- 
ing this characteristic size is best interpreted as the 
diameter of the Hll region (note that 6p oc porf^/t). 

3.1.1. Regime of Validity 

The above equations must be restricted to HlI re- 
gions that are still expanding and bounded within 
their GMC at the end of their ionizing lifetimes. If 
its luminosity is too low, the region will decelerate 
to CefT and stall before its driving stars burn out; as 
WM97 argue, such associations are too small to mat- 
ter (but see §|8|). If its luminosity is too high, the 
HlI region will envelop its entire cloud (rn > Rd), 
deforming the GMC into a cometary configuration 
(Bertoldi & McKee 1990) rocketing away from the 
association ( port Sz Spitzer 1955| ) . Since this reduces 
the rate of photoevaporation considerably, WM97 ar- 
gued that an appropriate maximum value for JM^est 
should be roughly the value predicted for a blister re- 
gion at twice the time required for its size to match 
the cloud radius. Using Whitworth (1979| )'s theory 
for the size scale, WM97 identified a maximum value 



5M, 



dest,max 



4.6 X 10"^ 



Nh. 



22 



1.5 



-3/8 



(20) 

Because the blister H II region radius derived in equa- 



tion (17) is smaller by a factor 1.6 than in Whitworth 



(T97E )'s theory, this upper limit corresponds to eval- 
at a time when rn 



uating equation ( jl^ ) at a time when rn = 0.94i?ci- 
Since this upper limit is quite uncertain, and since 
there is an ambiguity between the radius and diame- 
ter of the Hit region in WM97's argument, we shall 



simply adopt their value as given in equation (20). 
Correspondingly, we shall take an upper limit for the 
effective injection of momentum to be 



2cu6Md. 



est, max* 



(21) 



The axisymmetric calculations of Yorke et al. (198"9[ ) 
give a result that is only 6% lower, once we account 
for the differences between their ionized sound speed 
and recombination coefficient and ours (without these 
corrections, their result would be 17% higher). This 

^The solution for embedded regions was presented in 1995 by C. 



Note that these upper limits on dM^est and 5p are 
important for Hll regions if 



5'49 k, 63 



Nh. 



22 



1.5 



-3/8 



M, 



5/2 



cl.6 



3.7 Myr 



(22) 



F. McKee in lectures for Ay216 at U.C. Berkeley. 
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A more serious limitation on the theory presented 
here arises from the fact that we have ignored grav- 
ity in the evolution of an Hll region. If the GMC is 



through stellar ejecta rather than through their Hll 
regions; this is significant for the massive clusters 
forming in starburst environments ( [Tan fc McKee 



too d snse [tih > 140 cm '^), its free-fall time will be 200C ) and for the most luminous of Milky Way cluS' 



shorter than the typical ionizing lifetime of 3.7 Myr. 
This is typically true of clouds unless 



Mci > 3.7 X 10^ Mq. 



(23) 



For clouds below this limit, the orbital motion of an 
association within or about the cloud is likely to al- 
ter its Hll region, possibly by converting it into a 
cometary H II region of the type discussed by IRasi 



wala 1(19651 ), |Raga (1986| ) and |Raga et al. (199?! ). 



Lastly, we have not attempted to account for the 
inhomogeneities of molecular cloud material in the 
evolution of an Hll region (e.g., pyson et al. 1995| ), 
nor to instabilities that may develop during its ex- 
pansion ( Parcia-Segura fc Franco 1996| ). The finite 
porosity of the GMC (i.e., the interaction between 
Hll regions; see WM97) will be accounted for in an 
approximate manner in §^. 

3.2. Contribution from Stellar Winds and 
Supernovae 

We now consider the effects of stellar winds and 
supernovae on the evolution of an Hll region. Re- 
sults in this section are based on the stellar evolution 
code Starburst99 ( |Leitherer et al. 1999| ) , for which we 
have used a Scalo stellar initial mass function nor- 
malized as in MW97. In the first 3.7 million years 
after stars form, they inject roughly 78M^ km s""*^ of 
momentum: 40i?!>^M^km s~^ from protostellar winds, 
30M*km s~^ from main-sequence and evolved stars, 
and 8M^km s~^ in supernova ejecta. In the formation 
of a sufficiently massive stellar cluster, protostellar 
winds are mostly stopped within the self-gravitating 
clump from which the cluster arose ( Matzner fc Mc-| 
Kee 2000D and will affe ct this clump's dynamics 
( [Bertoldi fc McKee 19961 pVIatzner fc McKee 1999bD 
more than those of the surrounding cloud. The re- 
maining contribution is 



6pci ~ 38 



t 



3.7 Myr 



M^km s 



-1 



(24) 



If, in the first 3.7 Myr, this is smaller than the mo- 
mentum imparted by the HlI region, the ram pressure 
from stellar winds is lower than piiCjj. This means 
that the winds will be confined within the Hll region, 
while the ionizing stars shine, unless 



549 ^ (270, 210) X 



22 



1.5 



-1/2 



M 



1/6 
cl,6 



(25) 



for (blister, embedded) HlI regions, respectively. 
Sufficiently luminous associations couple primarily 



ters. 

Equation (^5|) accounts only for the momentum in 
material flung away from a star. More radial mo- 
mentum can potentially be generated in a pressurized 
bubble or blastwave that entrains ambient mass, as 
the radial momentum varies with energy E and total 
mass M as (5p ~ {EMY^"^ . Whereas in a thin shell 6p 
is conserved while E is not, in a fully adiabatic bub- 
ble or blastwave E is conserved while 5p increases as 
M does. Intermediate cases, such as pressure-driven 
snowplows ( Ostriker fc McKee 198^) a nd partially ra- 
diative bubbles ( Koo fc McKee 1992| ), involve a loss 
of energy but a gain of radial momentum. 

McKee et ai. (1984) have argued that, in the con- 
text of an inhomogeneous medium, stellar wind bub- 
bles are either confined within their HlI regions or 
made radiative by mass input from photoevaporat- 
ing clumps. In a blister HlI region, furthermore, hot 
gas can escape from the cloud. For these reasons we 
shall assume that stellar wind bubbles do not gen- 
erate momentum significantly in excess of the wind 
momentum itself. 

Supernovae merit special attention, as supernova 
remnants typically experience an adiabatic Sedov- 
Taylor phase {5p tx r^/^) followed by a pressure- 
driven snowplow phase {5p oc r^/^). For an up- 
per IMF cutoff of 120 Mq, the first supernovae 
explode 3.6 Myr after the onset of star forma- 
tion and thereafter occur every (3/S'49) Myr for 
the first million years, slowing to once per (5/S49) 
Myr thereafter. This frequency should be com- 
pared with the sound-crossing time of the HlI region, 
2(Mci,65'49)^/^(iVj/,22/1.5)"^/'^ Myr, and with the e- 
folding decay time of the ionizing luminosity, roughly 
2.5 Myr. In a rich association the delay between su- 
pernovae is the shortest of these, implying that su- 
pernovae blend together with the stellar winds. The 
above arguments indicate that the combined effects of 
supernovae and winds are negligible unless equation 
(P^) is satisfied. 

As a check, we have calculated the dynamics of the 
very first supernova remnant inside the HlI region, 
using the theory of Cioffi et al. (198'^ ). Because the 
progenitor mass is very high (~ 90 Mq in the presu- 
pernova state), the remnant becomes radiative before 
a comparable mass has been swept up, skipping the 



Sedov- Taylor phase (see also [Wheeler et al. 19"80D . 
The high progenitor mass also allows for a transition 
to a momentum-conserving snowplow phase while the 
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remnant is still expanding supersonically. The rem- our estimate Ces — 0.57urms is based on McKee 



nant becomes subsonic and merges with the Hll re- (1999| )'s formula for the mean pressure within GMCs. 



gion before striking its periphery if 



^49 > (2.0,2.7) X 



22 



V 1.5 



0.66 



(26) 



for (blister, embedded) regions. If this happens, the 
remnant adds little or no momentum to the Hll re- 
gion. This corroborates the above conclusion that 
supernovae are not significant for the momentum of 
the Hll region, except for regions smaller than the 
limit in equation (|2^) or larger than the limit in equa- 
tion (|25|). The possibility remains that supernovae 
in small associations contribute non-negligibly to the 
total momentum input and hence reduce the ioniz- 
ing flux and star formation rate we derive below by 
neglecting supernovae; this merits further study. 

3.3. Direct Disruption of Small Clouds? 
If a cluster were to deliver an impulse in excess of 



22 



ii4_ 



M^/q Mq km s-i 
(27) 

to its parent cloud in a time short compared to 
the cloud's dynamical time, the cloud would be dy- 
namically unbound. This is an unattainably large 
value for giant clouds with masses of ~ 10^ Mq, 
but smaller clouds might be disrupted in this man- 
ner. The Taurus- Auriga and Ophiucus clouds have 
Mcifi — 0.01, and may be susceptible to disruption. 
Unfortunately the likelihood and frequency of this 
process is difficult to assess, as it depends on several 
uncertain elements: (^Pmax 

(eq. |l|), the maximum 
size of a cluster that can form within a given cloud 
(WM97), and the importance of gravity for HlI re- 
gions in small clouds (§ |3.1.1| ). 

3.4. Turbulent Forcing Scale 

For turbulence driven by expanding shells, the rel- 
evant forcing scale is the radius of a shell once it has 
decelerated and become subsonic relative to Cefr- This 
is 



mm ■ 



3 5p 



(2, 4) X vrpoCcfT 



1/3 



cl 



(28) 



for (blister, embedded) regions, respectively. We do 
not allow rm to exceed the cloud radius; ignoring stel- 
lar winds, this limit becomes important when 



549 > (54, 220) X 



Ccflf 



0.57v, 



7/4 /^T 



^H,22 

1.5 



0.81 



M, 



2.1. 

cl,6' 

(29) 



How is the eff'ective forcing wavelength Ain related 
to the merging radius in equation (^)? Diametri- 
cally opposed regions of a disturbance move in oppo- 
site directions, and should be considered half a wave- 



length apart: this suggests Ai 



4rm. A similar 



conclusion follows from the fact that a shell's expan- 
sion could be considered a quarter cycle of oscillation 
(and its collapse, a second quarter cycle); for these 
reasons we set 



Ai, 



(30) 



and consider (pm an uncertain parameter of order 
unity. We shall find in §§5.1 and ^ that this is roughly 
consistent with the Galactic ionizing luminosity and 
star formation rate. Conversely, if future numerical 
simulations indicate a value of the product (j)m4' much 
less than unity, then feedback from H II regions would 
not explain the ionizing luminosity and star forma- 
tion rate in the inner Galaxy. This question should 



soon be addressed by simulations like those of Mac 
Low (2000|) and |Mac Low et al. (200l|) . Note that 
this requires Ceg and to be identified in the simu- 
lation volume, so that loss of energy by decelerating 
shells can be discriminated from turbulent dissipa- 
tion: combining equations (0), (|3|), (0), and (BH), 



3 

rms 



1 / 0.78MciW 

2 I 0.83^jinrmCefr 



(31) 



where refers to the creation of momentum by ex- 
plosive events in a simulation driven by point sources, 
and the coefficients 0.78 and 0.83 should be ad- 
justed to agree with the dissipation rate and degree 
of equipartition, respectively, in simulations of homo- 
geneously driven turbulence (as in eqs. Q and [^) 
that share global parameters like the ratio of mag- 
netic to gas pressure. The numerical determination 
of the product (/>m,'^ is essential to evaluate the theory 
presented here, or indeed to estimate stellar feedback 
in general. 

4. FEEDBACK FROM A POPULATION OF OB 
ASSOCIATIONS 

We now wish to account for the dynamical feed- 
back from the entire population of HlI regions that 
will exist within a given GMC. A similar project was 
undertaken by MW97, who modeled the luminosity 
function of HlI regions in the Galaxy, and WM97, 
who calculated the lifetimes of GMCs after hypoth- 
esizing how the Galactic population of HlI regions 
was distributed amongst GMCs. Our discussion fol- 
lows these prior works to the greatest degree possible, 
except that we shall solve self-consistently for the rate 
of star formation within a given GMC. 
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MW97 adopted a Scalo (1986 ) stellar initial mass 
function (IMF), with an upper cutoff at 120 Mq, nor- 
malized for a mean stellar mass (m^) = 0.51 Mq. 
They used the results of Vacca et al. (1996| ) for the 
ionizing fluxes and lifetimes of massive stars: aver- 
aged over the IMF, this gives a mean ionizing flux 
(549) = 8.9 X 10~^ and a mean ionizing lifetime 
{tras) = 3.7 Myr. They assumed that the fraction 
of OB associations born with more than A/"* stars, 
which we denote Fa{> A4), satisfies 



dFa 



1 



oc 



(32) 



from a lower limit of 100 up to a maximum of 5.5 x 10^ 
stars. With these limits, MW97 were able to fit sur- 
veys of luminous Hll regions (e.g., Kennicutt et a\\ 
1989|]|, the total star formation rate in the Galaxy, 
the birthrate of nearby associations, and the Galac- 
tic recombination radiation. 

A cluster's ionizing luminosity is dominated by 
its most massive members; this implies a distinc- 
tion between rich clusters (^49 ^ 10; Kennicutt 



et al. 1989), which are populous enough to sample 
the IMF up to the upper cutoff, and poor clusters, 
which are not. The ionizing luminosity of a poor 
cluster reflects its most massive member, and rises 
only statistically with A/"*. The cluster's ionizing life- 
time likewise reflects this star's main-sequence life- 
time (tms oc 5^9°-^^; [Vacca et al. 1996| , MW97). In 
contrast, a rich cluster's ionizing luminosity and life- 
time represent averages over the IMF: 5*49 = (549) A/* 
and ti = (tins) 5 respectively. This distinction intro- 
duces a turnover in the cluster luminosity function 
Fa{> Sig) at the boundary between rich and poor 
clusters: dFa/dln S^g oc S^g for rich clusters, but 
statistical fluctuations cause a flattening of -Fa(> ^49) 
for poor clusters (iKennicutt et al. 1989| and MW97). 
Associations near this turnover dominate the feed- 
back from HlI regions on GMCs, because 5p (and 

for blister regions, (JM^est) scale as S^^t^J'^. For in- 
stance, WM97 found that half of all photoevaporation 
in clouds with Mf.\fi > 0.1 is accomplished by regions 
with 549 < 3.7. 

It will be useful to define the cluster-weighted mean 
of a quantity x as 



(33) 



For MW97's Galactic Hll region luminosity function 
the mean mass per association is {M^)^ = 440 M©, 
corresponding to (549) = 0.77. Other useful aver- 
ages are listed in table ^. 

The Galactic population of GMCs, Mc\{> M^), 
satisfies d/Vd/ciln Md oc Mj" where a ~ 0.6 from 



an undetermined lower limit to an upper limit of 
6 X 10*^ Mq (MW97). MW97 argue that a given cloud 
cannot make arbitrarily large OB associations. Tak- 
ing the maximum cluster mass to be 10% of the cloud 
mass, they derive the upper limit 

(34) 



^49 < 5«,49(Mci) = min(490, 172Md,6) 



which we also adopt. Hll regions must therefore oc- 
cur within a given GMC in different proportions than 
they are found in the galaxy. WM97 construct a 
cloud's luminosity function -Fa,A/(> 5*49) by assum- 
ing: 1. no clusters form above S'u,49(Mci); 2. since 
GMCs give birth to HlI regions, the Galactic lumi- 
nosity function must equal the sum of all GMCs' lu- 
minosity functions; and 3. within the clouds that can 
form a given size of cluster, the birthrate of those 
clusters is proportional to cloud mass. These assump- 
tions led to MW97's equation (23) for the population 
of associations forming within a given cloud. Below, 
we shall solve for the star formation rate in a given 
cloud under the hypothesis that it is supported by Hll 
regions forming within it; this requires that assump- 
tion (3) be dropped. If the star formation rate within 
clouds scales as oc M^, a derivation analogous to 
MW97's equation (16) gives 



Fa,Mi> Of 



H[K,uiM,i) - 



1 - (AA^/1.2 X 106)/3 



(35) 

where H{x) = (1,0) for (x > 0,x < 0) is the step 
function and A/;,«(Mci) = SuM^d)/ {sAg)- MW97 
assumed (3=1; below we derive /3 = 1.32 for the mass 
range of interest, and in §0 we find (3 = 1.38 after ac- 
counting for the interaction between HlI regions. The 
above birthrate distribution is used in Monte Carlo 
simulations (as described by MW97) to produce the 
figures in this section. However, the following argu- 
ment indicates that the difference between Fa and 
Fa^M is unimportant for clouds in the mass range 
where our theory is valid. 

4.1. The Giant-Cloud Approximation 

For analytical estimates it is useful to note that the 
associations responsible for most of the mass ejection 
and energy injection (S^g ~ 3.7) are much smaller, for 
giant clouds, than those associations whose existence 
or behavior are affected by their sizes or the finite 
sizes of their parent clouds. Specifically, we noted 
m §|T| that our neglect of cloud gravity is only 
valid for clouds with M^fi > 0.37. For such clouds, 
Su 49 > 64; the HlI regions whose maximum sizes (eq. 
[2^ ) and merging radii (eq. [^) approach Rd have 
5*49 > 5, increasing rapidly with Mci; and those for 
which winds and supernovae are important (eq. ||25|| ) 
have 549 > 180. These facts have two implications: 
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Table 1 

Moments of the Galactic luminosity function of Hi: regions 



Variable 



J 



5^9(^^/3.7 Myr)' 



Eqn. 



{6p)^ , (<5Afdcst) 

tdO 

Qo 



4/7 
16/21 
1/14 

3/7 



9/7 
12/7 
9/7 
19/7 



0.499 
0.572 
1.856 
2.128 



(45) 
(46) 



m 
m 



Note. — Averages over the Galactic initial luminosity function of OB associations, derived from a Scale 
(1986) IMF and a Monte-Carlo simulation as described by MW97. 



1. A cloud's luminosity function Fa^Mi> S^g) only 
differs from the Galactic form Fa{> S4g) for luminosi- 
ties near 5^4^49 (Md), so we may neglect this difference 
in analytical estimates; and 2. We may assume that 
Hll regions are entirely contained within their par- 
ent clouds. These approximations, to which we refer 
collectively as the giant- cloud approximation, are the 
same that led to MW97's equation (40). 

For clouds with M^fi < 3.7 it might be more appro- 
priate to take a small-cloud approximation in which 
all Hit regions are assumed to outgrow their clouds: 



SMdest = 5Mdest,max, Sp = 6p„ 



and 



Rr 



(cf. MW97 eq. [41]). This would be very uncertain, 
however, both because our estimates of 6M([est max 
and 5pmax are uncertain and also because we have 
neglected cloud gravity, which should be important 
when Mci < 3.7 x 10^ Mq. 

The performance of the giant-cloud approximation 
is illustrated in figure [l|, where we plot its prediction 
for {5p)^ [dashed line, eq. ||T^) against the results 
of a Monte-Carlo simulation (solid line). The simu- 
lation accounts for the difference between the lumi- 
nosity distribution of OB associations within a given 
CMC (Fa^M, eq. ||35||) and that within the Galaxy 
as a whole {Fa, eq. p2| ), whereas this distinction is 
neglected in the giant-cloud approximation. This is 
the p rimary cause for the difference between the two 



A comparison between {5p)^ and {M^)^ gives the ef- 
fective input of momentum per stellar mass averaged 
over Hit regions. In the giant-cloud approximation, 
an integral of equation (O) over Fa{> S^q) gives 



260 

\ 1.5 



-3/14 



mKI" km s-i X (M.)(^6) 



The characteristic velocity (momentum per stellar 
mass) identified in this equation is much larger than 
the corresponding coefficient 5pw/nn-^ = 40i;^>^km s~^ 
that we have estimated for protostellar winds in equa- 
tion (|5|), or the value 38 km s~^ estimated for main- 
sequence and evolved star ejecta in equation (24). 
Since equation (p6|) applies to both blister-type and 



embedded Hll regions, it implies that HlI regions are 
more important than the combined effects of proto- 
stellar winds, stellar winds, and supernovae in driving 
turbulence within GMCs, regardless of the efficiency 
of star formation. Figure || makes this point graphi- 
cally. 

Our estimate of (^Pdest,max indicates that HlI re- 
gions dominate over protostellar winds so long as 
Mci ^ 4x 10^ Mq. This, in turn, suggests that proto- 
stellar winds may support small clouds and the self- 
gravitating clumps within GMCs ( [Bertoldi fc McKee 



curve b for Mh > 3.7 x 10^ MfT,, to the right of the 1996| ; [Matzner fc McKee 1999bD , whereas HlI regions 



vertical dotted line; for instance, the kink in the sim- 
ulation at Mc\fi = 3 is caused by the kink in Fa^M 
in equation (^). The simulation also accounts for 
finite-size effects (<5pmax) eq. [^); these dominate 
the difference between the curves to the left of the 
dotted line. However, our neglect of cloud gravity is 
not valid in that region (eq. |l23|| ). 

4.2. Momentum Generation 



support the GMCs themselves. 



4.3. Mass Ejection 



If we assume that all Hll regions evolve into blis- 
ter regions, as have Whitworth (1979| ), [Franco et al 
19941 ), and WM97, we arrive at an upper limit to the 



amount of mass than can be removed from a CMC 
by photoevaporation. Since 5p = 2cii(5Mdest for a 
blister region, equation (|3^) gives, in the giant-cloud 
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Fig. 1. — Cloud momentum generated by HII regions per stellar mass formed, and an illustration of the giant-cloud approxima- 
tion. The theory for HII regions presented here does not account for gravity and is therefore only valid for clouds to the right of 
the vertical dotted hue ~ those whose free-fall times are longer than the ionizing lifetimes of OB associations (eq. [^). The solid 
line is the result of Monte-Carlo simulations that account for the difference between the Galactic H II region luminosity function 
and what is expected within a given cloud (eq. [^; see also MW97). These simulations also account for the effect of finite cloud 
size. The dashed line represents the giant-cloud approximation in which both the corrections to the Galactic luminosity function 
and the finite-size effects are ignored. To the right of the dotted line, differences between the two curves result primarily from the 
luminosity function (and are relatively minor): for instance, the kink at M^i = 3 x 10® Mq divides clouds that can harbor the 
largest Galactic OB associations from those that cannot. Finite cloud-size effects dominate to the left of the dotted line, where our 
model of HII regions is already invalidated by gravity. The dash-dot lines show the momenta of protostellar and main-sequence 
winds in the first 3.7 Myr of a cluster's life. 
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approximation, 



{SM, 



dest/ 



14 



22 



V 1.5 



-3/14 



Mill' X (M.)„ , (37) 



implying a typical instantaneous star formation effi- 
ciency 



~ 6.7% X 



3/14 



-1/14 



cl,6 



(<5Mdest)„ + (M,)„ V 1.5 

(38) 

SFEtot and e are related by SFEtot = (s^^) 
( [Matzner fc McKee 200Cl| ); however, a sudden disrup- 
tion of the cloud (a precipitous drop of e, as discussed 
in §3^) would reduce SFEtot below the typical value 
of £. 

The lines marked A in figure Q compare equation 
(^) with the results of a Monte-Carlo simulation 
which accounts for the birthrates of OB associations 



within individual clouds as given by equation (35). 
The lines marked B account self-consistently for the 
frequent interaction between Hll regions, in a manner 
we discuss in §0. 

Unfortunately, there are insufficient observations of 
clouds in the mass range 3.7 X 10^ Mq < Mci < 
6 X 10® Mq to test equation (^) directly. Observa- 
tions of nearby clouds, which are too small for an ap- 
plication of our theory, indicate a significantly lower 
star formation efficiency than an extrapolation of our 
theory would predict: for instance, observations of 
nearby small clouds ( [Evans &: Lada 1991 ) indicate 
that only ~ 1% of the mass is stellar and only ~ 3% 
will ever be. Besides gravity, a number of other effects 
may be important for clouds this small; WM97 have 
emphasized their disruption by Hll regions that out- 



grow their boundaries (see also Elmegreen 1979), as 
well as the likelihood that this disruption will render 
them susceptible to photodissociation. 

5. IONIZING LUMINOSITY AND STAR FORMATION 

RATE 

As a cloud's turbulence decays it must be replen- 
ished, if virial balance is to be maintained. Energy 
can be derived from gravitational contraction, winds 
and supernovae, and Hll regions. Contraction is 
problematic as it would imply a collapse rate close 
to free-fall, hence rapid star formation; this may hold 
for the formation of stellar associations, but cannot 



for entire GMCs (Zuckerman fc Palmer 1974) 



We have shown that Hll regions are the most ef- 
fective of the remaining sources. Accordingly we will 
match the loss of turbulent momentum with its regen- 
eration in Hll regions. We must account for a pop- 
ulation of H II regions with different luminosities and 
therefore different values of 6p, r^, and Ain. Equa- 
tions (|^) and (^) predict that energy decays in a time 



proportional to Ain. It is most consistent to assume 
that the contribution from each type of Hll region 
decays independently on its own timescale; the to- 
tal turbulent energy is then the sum of these decay- 
ing contributions. Taking the decay time from equa- 
tion (^), and taking the energy input SEiy^j.]^{6p) from 
equation (^, 

EtuTh = X! "^-^turb X (decay time) x (rate) 

•S'49 

'^'^^■^'^ (Ain(M5p>. , (39) 



2AVr 



and 



where Xmi^p) is given by equations (p^) and (| 
Ma is the formation rate of associations. 

Applied to a CMC, equation (^) must be consis- 
tent with the observed kinetic energy of the cloud and 
the expected degree of equipartition between kinetic 
and magnetic energy (eq. 0). Equating expressions 
(0) and (p9|) gives the formation rate of associations: 



1.9MciM 



3 

rms 



(40) 



If one were to vary the ionizing flux of all stars in 
the IMF by the same factor, then Ma would vary as 

l/(Ain5p) oc i5p)-'/^ cx (549):''/'' oc (549)"''/''. 
This scaling holds also for the star formation rate 
= {Mi,)^Ma- But, the mean ionizing flux pro- 
duced by these associations would vary much less: 

S49,t{Mci) = {SA9)a (*)ms^a OC (549)^/''. This is not 

surprising, as the ionizing photons are directly re- 
sponsible for sustaining equilibrium. 

A cloud's ionizing luminosity is therefore quantity 
most tightly constrained by the assumption of equi- 
librium: in the giant-cloud approximation. 



(37,53) fNn, 22 



<t>m{(^/l.6) V 1.5 



1.37 



CefT 



0.57Vr 



XM, 



37/28 



cl,6 



-2/3 



(41) 



for Hll regions that are primarily in the (embedded, 
blister) state, respectively. This equation will be ad- 
justed for the interaction between Hll regions in §^. 
To compare with observations of individual clouds 
one must account for the leakage of ionizing photons; 
MW97 estimate that ~ 70% escape the observed Hll 
regions. 

The star formation rate is = {M^,)^Ma- Using 
the giant-cloud approximation and MW97's IMF, 
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Fig. 2. — Star formation efficiency in giant clouds, assuming photoevaporation is the dominant disruption mechanism. The 
curves marked A are not corrected for the interaction between If II regions, whereas the curves marked B account self-consistently 
for finite cloud porosity (§^. The solid lines are Monte-Carlo calculations based on MW97's model for the luminosity functions 
of OB associations forming within a given cloud; effects due to finit e cloud size (eqs. [pO|, fpl[ , and [Q) are included. The dashed 
lines are analytical results of the giant-cloud approximation (§ [l.l[ ): equations (Bq) anc (Isfl) for cases A and B, respectively. In 
this and subsequent figures we assume that all HII regions enter the blister stage (see 



for (embedded, blister) regions. Recall that oc 
-16/21 -g predicted less robustly than 
oc (,349)^^^^ due to uncertainties in the 



mean stellar mass and ionizing luminosity. [Kennicutt 
et al l (1994D find that the ratio (m*) /((S49) (tms)) 
varies by about 50% among the forms of the IMF 
and sets of stellar tracks they consider; the predicted 
is therefore uncertfiin by at least this amnimt 



Thfe mean ionizing luminosity and lifetime de- 
pend on the metallicity Z of the stellar population, 
and these effects cause a shift in the star forma- 
tion rate relative to equation (^2|) even if the IMF 
remains constant. Numerical integrations of equa- 



tion (15), performed for Starburst99 populations of 
various Z, exhibit {Sp)^^ / {M^)^ oc jniplying 
Mv, oc Z°-i2j\rL37_ g^j^y situation where M89's the- 
ory holds, Nh oc Z^^ and therefore oc Z^^''^^. 

Our hypothesis that the star formation rate is de- 
termined by a balance between the dissipation of tur- 
bulence and its regeneration in Hll regions can only 
make sense if it predicts more than one Hll region 
per dynamical time of the cloud. This is true for all 
the clouds within the regime of validity of our theory, 

as we find that Mats ^ (13, 19) x Nl^^!^f {1.6/[(f)^m(p]) 
so long as M^fi > 0.37. 
In figure |3| we plot tg^: = Md/Mi,, the time scale 



for gas to be converted into stars. Equation (|42|) is 



shown as the dashed curve marked A. The solid curve 
marked A is the results of a Monte-Carlo simulation 
which incorporates equation ( |32|) and the effects of 
finite cloud size (eqs. |^] and ||2^). The curves 
marked B are adjusted to account for the interaction 
between Hll regions, as described in §0. Also plotted 
{dash- dot line) is value of tg* appropriate t o the entire 
inner Galaxy, and the range suggested by Carpentei 
(200C| )'s observations of clouds around 2 x 10^ Mq, 



190-440 Myr (corresponding to the acceptable age 
range 3-7 Myr for the embedded population; see be- 
low). 

5.1. Galactic Ionizing Luminosity and Star 
Formation Rate 

Equations ( pl) and (^) refer only to a single molec- 
ular cloud. Integrating over the GMC population 
for the inner Galaxy, we arrive at the total inner- 
Galactic ionizing luminosity St and star formation 
rate M^t- WM97 model the GMC population of the 
inner galaxy as dNd/dlnM^ ~ 63(6/McL6)"°-^; set- 
ting Nh,22 = 1-5 as observed, and considering only 
giant clouds with M^ifi > 3.7 for which we are rela- 
tively confident of the star formation rate, we find 
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Fig. 3. — The time scale for molecular gas to be converted into stars, i^* = A'/ci/Af*. Here, as in the other figures, the vertical 
dotted line delimits the lower mass limit of validity of our theory (eg. Pq]); the solid lines are Monte-Carlo calculations. The 
dashed lines are the giant-cloud approximation, equations (^2|) and ( |4s| ) for A and B, respectively. The curves marked A neglect 
interactions between H II regions, whereas those marked B account self-consistently for these interactions in the manner derived in 
The value of tg* for the entire Gala xv is plotted a s the dash-dot line; this corresponds to a to tal molecular ma ss of 10^ Mq and 
star formation rate of 5.1 Mq per year (Mezger 1987). Also plotted is the value of tg* implied by Carpenter (200C)'s determination 
that the nearby Perseus, Orion A, Orion B, and Mon R2 clouds (10* Mq < M^i < 3 x lO" Mq) have converted ~ 1.6% of their 
mass into stars in the past 3-7 Myr. In this figure the product 0™,'/' is taken to be 1.6 x 1.2. 
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These results would be 15% higher if we were to ex- 
trapolate the giant-cloud approximation below the 
lower limit of Mc\fi = 0.37. As before, St is much 
better constrained than M^,t by our hypothesis that 
Hll regions dominate feedback. 

We may use observations to check that clouds with 
Mci^G < 0.37, whose evolution we cannot address, do 
not dominate the Galactic star formation rate or ion- 
izing flux. Consider Garpenter (2000| )'s observations 
of the Perseus, Orion A, Orion B, and Mon R2 clouds 
(Mci ~ 1 — 3x 10'^ Mq for each). Parpenter estimated 
the completeness of the 2MASS survey and the result- 
ing stellar fraction in these clouds, a slowly increasing 
function of the assumed age of the observed stellar 



population. Carpenter argues that that these clouds 
have created a stellar mass equal to ~ 1.6% of their 
own mass in a period between 3 and 7 Myr.'^ Extrap- 
olating this star formation rate per unit molecular 
mass to ah clouds with Mdfi < 0.37, using MW97's 
cloud mass function, gives a range of 0.7 to 1.7 Mq 
per year. This estimate should be considered very 
uncertain, as it assumes oc for small clouds; 
however it indicates that clouds in the mass range 
Mci < 3.7 X 10^ Mq, which constitute one third of 
the molecular mass, most likely do not dominate the 
Galactic star formation rate. 

MW97 estimate the ionizing flux of the inner 
Galaxy at 2.6 x 10^'^ s~^ (assuming all ionizing pho- 
tons are caught within the disk ) and derive a star 
formation rate of 4.0 Mq 'yv~^. [ Mezger (1987 ) esti- 
mates the star formation rate in the inner Galaxy to 
be 5.1 Mq yr~^.] Equations (]43| ) and (|4^ ) are com- 
parable to this value, with i?i> ~ 1.6 as suggested by 
McKee (1999 ), and with (j)^ — 1-1 (close to unity 
as suggested in §|3.4| ). Given the approximations that 
led to equation ([44D , this degree of agreement is quite 
remarkable; it supports our proposition that feedback 
from Hll regions regulates the rate of star formation 
in giant molecular clouds. 

6. LIFETIMES OF GMCS 

We may combine the star formation rate given in 
equation ( p2[ ) with the upper limit for mass evapora- 
tion from GMCs given in equation (|37|), to arrive at a 
lower limit for the lifetime of GMCs assuming all H II 
regions evolve into blister regions. In the giant-cloud 
approximation, this gives 
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for Mci > 3.7 X 10^ Mq. A significant population of 
Hll regions that remain embedded rather than evolv- 
ing into blister regions would extend this lifetime, but 
we argue in ^ that such a population is not likely. 
Another effect that extends cloud lifetimes relative 
to trfo is the finite porosity of Hll regions within their 
volumes (MW97), for which we account in §|^. 
Note that the estimate of the cloud destruc- 



roughly 2.6M 



-4/7 



cl,6 



tion time in equation (^) is 
times longer than the cloud's free-fall time; the 
two time scales are equal for a mass M^i = 5.3 x 
106(iV//,22/1.5)-13/12 Mq, very close to the upper 
mass hmit for Galactic GMCs (~ 6 x 10^ Mq; 
WM97). This raises the possibility that the upper 
mass limit derives from the difficulty of assembling 
an object that destroys itself rapidly compared to its 
free-fall time scale, a hypothesis that should be tested 
against extragalactic observations. We refine this ar- 
gument in §1^, where we adjust the destruction time 
for finite cloud porosity. (But note, a possible coun- 
terargument is raised in §^.) 

Cloud destruction timescales are plotted in figure 
1^, where our Monte-Carlo simulation [solid lines) is 
compared with the results of the giant-cloud approx- 
imation as given by equations (45) and (^) for cases 
A and B, respectively. The curves marked B are ad- 
justed for the interaction between Hll regions (^. 
The free-fall time is also plotted, making visible the 
crisis of rapid destruction for clouds above the Galac- 
tic upper mass limit. 

Our equation ( p5|) is based on the same physics 
that led MW97 to estimate a destruction time of 20- 
25 Myr for clouds with M^fi > 0.1. Our estimate is 
shorter because MW97 assumed oc Md, whereas 
an energetic balance requires oc M^^^ in our giant- 
cloud approximation; this redistribution shortens the 
lives of the massive clouds to which our theory ap- 
plies. 

7. EFFECT OF FINITE POROSITY 

The frequency with which Hll regions form within 
GMCs implies that they are likely to interact 
(WM97). Defining the porosity Q as the time- 
averaged volume filling factor of Hll regions, Q = 
AfaJ J^'''^[Msh{t)/Mci]dtdFa, we estimate from the 
above equations 
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Fig. 4. — Destruction time scales tdo (case A, ignoring cloud porosity) and td (case B, accounting for finite porosity) for giant 
clouds. The solid lines are Monte-Carlo calculations; the dash-dot lines, which result from the giant-cloud approximation, repre- 
sent equations ( ^sj ) and (|^ for cases A and B, respectively. Also plotted is the clouds' free-fall time fg, which approaches their 
destruction time at the upper mass limit for Galactic GMCs. 
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Fig. 5. — The porosity of HIT regions within GMCs. Plotted are the initial estimate {Qo, case ^4), which does not account 
self-consistently for the effect of porosity on star formation rate, and the revised estimate {Q, case B), which does. The dashed 
lines result from the giant-cloud approximation; case A is given by equation (^^. 



17 



for (blister, embedded) regions respectively. This ex- 
pression is not yet fully self-consistent; we find Q in 
terms of Qq below. The large value of Qq for embed- 
ded regions implies that fully embedded Hll regions 
are ruled out, as Hll regions will percolate and al- 
low their gas to vent. This reinforces our neglect of 
embedded regions in equation (^). 

WM97 argue that the effect of interactions is equiv- 
alent to regrouping stars into non-interacting re- 
gions that are larger by the factor (1 -|- Q), i.e.. 
Mi = Ma/{1 + Q) and = (1 + Q)549. For 
each of the regrouped regions, 

6M^^^^ = (1 + Q)'^'''^'^Mdest, and, for the giant-cloud 
approximation (r^ < -Rcl), = (1 + Q)'^^'^^hn- 
Whereas WM97 assumed a star formation rate within 
a given cloud, our theory demands that we determine 
the star formation rate self-consistently by matching 
turbulent decay with driving by Hll regions. Apply- 
ing these transformations to equation (E^) and elim- 
inating Mi, we find that Ma = {I + Q)^^MaO where 
Mao is the uncorrected value in equation (|40|) . But, 
the porosity is proportional the star formation rate: 
Q = QoMa/MaO- A self-consistent value of Q there- 
fore satisfies 

Q(l + Q)-^/2i = Qo, (47) 

so that Q ^ Qo when Qo < Ij but Q ~ Qq^^^^ 
when Qo ^ 1. We find that the useful relation 
{1 + Q) = 1.85Qo holds within 1% in the cloud mass 
range of interest. The analytical (equation |^6|) and 
numerical values of Qo and Q are plotted in figure ^. 
To account self-consistently for finite porosity, we 
incroQoo Ma, 6'49,r(j\fd) and by (1 I Q)^/^^ 



must 

relati Vo to oquationo 



:|lTD, and 



mcroaoo 



the cloud lifetimes by (l + Q)^''^^ relative to equation 
(^), and increase the stellar mass per ejected mass 
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(P|). The revised equaLiuiis are 
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respectively. 



The corrected lifetime is 



1.2{4>m(p/lW-^ {Md,6/6)-^-^ts, i.e., just over one 
free-fall time at the upper mass limit for Galac- 
tic GMCs. The total ionizing luminosity and 
star formation rate in the inner Galaxy due to 
clouds in the mass range 3.7 X 10^ Mq < Mci < 
6x10*^ Mq become 4.5 X 10^3 (0m (A/1-6)"^^/^^ s"^ and 
7.0((/)m(/'/l-6)~^^/^^ Mq yr~^, respectively. These are 
consistent with the observed rate if i;^ ~ 1.6 and 
0m ^ 1-6. 

Thus, our estimates of the Galaxy's ionizing lumi- 
nosity and star formation rate (with (pm'P — 1 x 1.6) 
are comparable to the values adopted by MW97 if 
cloud porosity is not accounted for (consistency re- 
quiring (pm4' — 1-1 ^ 1-6) and somewhat higher (con- 
sistency requiring (prn4' — 1-6 x 1.6) after the porosity 
correction is applied. This level of agreement sug- 
gests H II regions are indeed responsible for maintain- 
ing energetic equilibrium within GMCs. What might 
explain the remaining discrepancy? First, numerical 
estimates of turbulent decay rates may decrease as 
studies improve in resolution and less diffusive codes 
are used, o r as turbulent anisotropy is included ( phoj 



et al. 2001). Second, the Galaxy's ionizing luminos- 



ity could be higher than MW97 found, if a signifi- 
cant fraction of ionizing photons escape the Galactic 
disk (unlikely, but controversial; Bland-Hawthorn & 
Putman 2001] ) - but conversely, additional ionization 



from sources like supernova remnants ( Slavin et al 
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-o.832 00q ) would increase the discrepancy. Third, the dis- 



crepancy results from our method of accounting for 
finite porosity, a crude approximation when Q > 1. 
(48)Fourth, we were not entirely able to exclude a contri- 
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bution of momentum from super novae in §3.2, espe- 
cially for small associations. 

Lastly, we have followed MW97 in approximating 
clusters' ionizing luminosity as a step function of du- 
ration (tms) = 3.7 Myr. In reality, a numerical inte- 
gration of equation ( |l5| ) using the ionization history 
predicted by a Starburst99 synthesis indicates that 
5p continues to grow (albeit more slowly), gaining 
another 50% after 6.6 Myr. This should effectively 
increase (tms) in the largest clouds, whose free-fall 
times are long enough to permit this expansion. All 
of these topics merit further study. 
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8. CONCLUSIONS 
The main results of this paper are as fohows. 



1. 



Hll regions are the most plentiful sources of 
energy for the turbulence within giant molecu- 
lar clouds; they are more significant than the 
combined effects of protostellar winds, main- 
sequence and evolved-star winds, and super- 
novae. This result was indicated by the low 
efficiency of star formation in GMCs in §2.1 and 
demonstrated on more general grounds in ^^7^ 



and 4.2 



2. The input of turbulent energy by Hll regions 
occurs on scales comparable to, but somewhat 
smaller than, the cloud radius. Large-scale forc- 
ing minimizes the rate of turbulent decay, as 
recently emphasized by Basu fc Murali (2001 ). 



3. A balance between turbulent decay and the re- 
generation of turbulence by Hll regions allows 
a prediction of the stellar ionizing luminosity 
and (less robustly) the star formation rate. We 
present these results for clouds in the mass 
range 3.7 x 10^ M© < Md < 6 x 10^ Mq, in 
which the stellar ionizing lifetime is briefer than 
the free-fall time. The results are roughly con- 
sistent with the total ionizing luminosity and 
star formation rate of the inner Galaxy, pro- 
vided that future numerical simulations verify 
our estimate of the coupling between momen- 
tum input and turbulent energy (eq. ) . Our 
estimate of the ionizing flux in the inner Galaxy 
is ~ 70% higher than the observed value; we list 
in ^ several possible resolutions of this discrep- 
ancy. 

4. Because Hll regions also evaporate their clouds, 
an energetic balance also implies a rate of pho- 
toevaporation and a destruction time scale for 
GMCs. We derive a destruction time of 17 to 
24 times (Mci/lO^Mo)-^^ million years, some- 
what shorter than was found by WM97. 

5. The upper mass limit for Milky Way GMCs 
most likely derives from the difficulty of assem- 
bling an object that destroys itself in a single 
crossing time. However, less massive clouds (at 
least down to 3.7 x 10^ -^0) survive for many 
crossing times, produce many HlI regions per 
crossing time, and are in both energetic and 
dynamical equilibrium. 

6. So long as there exists a minimum optical depth 
or column density required for star formation, 
the vigorous energetic feedback by Hll regions 
provides a mechanism for the maintenance of 



cloud column densities near the critical value, 
and hence, for the CMC line width-size rela- 
tion. Massive clouds therefore follow the sce- 
nario proposed by M89, but with Hll regions 
rather than protostellar winds as the primary 
agents of feedback. This conclusion is not as- 
sured for clouds less massive than 3.7 X 10^ Mq, 
for which cloud gravity must be considered in 
the dynamics of Hll regions. 

The theory we have presented is robust, in the sense 
that it applies to massive GMCs (containing the bulk 
of the molecular mass) and derives from the flat- 
tening of the luminosity function between rich and 
poor OB associations - a product solely of the up- 
per mass limit for stars - rather than the detailed 
birthrate distribution of OB associations. So long as 
a galaxy's molecular mass is concentrated in the most 
massive clouds, and so long as the birthrate of stel- 
lar associations drops with ionizing luminosity more 

steeply than d/Va/filn ^49 oc , energetic equilib- 
rium within molecular clouds determines its ionizing 
flux. This assertion can be tested by extragalactic 
observations. Starbursts may however be an excep- 
tion to this rule if, as discussed below, their GMCs 
are not in equilibrium. 

Other observational tests of the theory presented 
here include the variation of ionizing luminosity and 
star formation rate (figure ^ with cloud mass and 
the high porosity of HII regions (figure ^) for massive 
GMCs. The variations of these quantities with mean 
cloud column density can potentially be tested by ob- 
servations of the SMC. The star formation efficiency 
(figure ^) and cloud lifetime (figure ^) will be more 
difficult to verify. 

The inhomogeneity of giant molecular clouds is the 
greatest source of uncertainty in the present work. 
The interaction of Hll regions ~ one source of inho- 
mogeneity - was treated in an approximate manner 
in §^, but future work must treat the interaction of 
Hll regions with a realistic background cloud. 

We have noted that clouds more massive than 
6 X 10^ Mq should not form in the Milky Way be- 
cause they would be disrupted by Hll regions in 
the time needed to assemble them. However, we 
must also note that clouds more massive than about 
3.7 X 106(A^H,22/l-5)"^ Mq may be incapable of driv- 
ing champagne flows because their escape velocities 
exceed the exhaust velocity (^\\c\\ = 19.4 km s^^ of 
ionized gas. This fact does not affect our derivation 
of the star formation rate in such clouds, but it does 
call into question our derivation of the lifetime for 
the most massive clouds (and hence our suggestion 
for the origin of the upper mass limit). Further work 
will be needed to resolve this issue. 
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Similarly, clouds more massive than about 1.9 x number of starbursts; this is sufficient to crush all 

10'''(iVj^_22/l-5)~^ Mq cannot have supersonic Hll GMCs in the mass range > 3.7 x 10^ Mq. 
regions, as their effective sound speeds exceed 
10 km s~^, the sound speed of ionized gas. Objects 

created in this state or pushed into it by an increase It is a pleasure to thank Chris McKee, Jonathan 

in external pressure can neither be disrupted nor even Wilhams, and the referee, Mordecai-Mark Mac Low, 

supported by photoionization, and must collapse as for substantial comments that led to improvements of 

rapidly as their turbulence decays; this is a recipe for this paper. I am also grateful to Peter Goldreich, Rob 

the efficient production of massive star clusters ob- Kennicutt, Peter Martin, Eve Ostriker, Nick Scoville, 

served to occur in starburst galaxies. For instance, and Jim Stone for insightful discussions and sugges- 
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